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HYDROPHONE  SENSOR  HAVING  A CARDIOID  DIRECTIVITY  PATTERN 


Introduction 

To  meet  the  demands  for  specialized  underwater  acoustic  measure- 
ments by  the  Navy,  the  Underwater  Sound  Reference  Division  (USRD)  of  the 
Naval  Research  Laboratory  is  continually  evaluating  new  hydrophone  sen- 
sor configurations.  These  measurements  often  require  hydrophones  with 
unique  directivity  requirements. 

The  USRD  was  asked  to  develop  a hydrophone  having  a directivity 
pattern  within  a limited  frequency  range  that  is  omnidirectional  in  the 
space  below  its  horizontal  plane  and  with  reduced  sensitivity  in  the 
space  above  this  plane.  This  requirement  results  in  a hydrophone  that 
is  less  sensitive  to  the  interference  of  water  surface  reflections  and 
acoustic  noise  due  to  sea  state  conditions.  The  hydrophone  was  also 
required  to  be  omnidirectional  in  its  horizontal  plane. 

The  purpose  of  this  report  is  to  present  an  electroacoustic  anal- 
ysis and  comparative  performance  of  a hydrophone  sensor  configuration 
designed  to  fulfill  this  specific  directivity  need.  Also,  a design 
procedure  is  shown  to  configure  a sensor  of  this  type. 


Scope  And  Approach 

The  scope  of  this  report  is  based  on  a request  fcfr  the  USRD  to 
develop  a hydrophone  having  constant  sensitivity  at  any  angle  below 
the  horizontal  plane  and  being  optimized  to  have  minimum  sensitivity 
above  this  plane  in  the  frequency  range  of  14.5  to  19.5  kHz. 

The  experimental  evaluation  of  a known  hydrophone  sensor  config- 
uration having  this  directivity  was  completed.  Various  sensors  based 
on  this  configuration  were  also  experimentally  evaluated.  A concur- 
rent electroacoustic  analysis  was  initiated,  and  the  resulting  theo- 
retical directivity  was  compared  to  the  experimental.  This  was  done 
so  that  an  accurate  design  procedure  could  be  developed  for  this  sensor 
configuration  to  be  applied  to  developing  a sensor  fulfilling  the 
directivity  specifications. 


Results 

Sensor  Configuration 

The  sensor  element  configuration1  analyzed  is  shown  in  Figure  1. 
It  consists  of  an  axially-poled  piezoceramic  tube  with  a concentric 


*A.  C.  Tims  and  T.  A.  Henriquez,  Patent  No.  3,947,802,  March  30,  1976. 

Note:  Manuscript  submitted  July  13,  1977. 
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circular  piston  cemented  to  each  end.  This  assembly  is  sealed  concen- 
trically by  0-rings  inside  a precision-bored,  air-filled,  tubular 
housing  made  of  aluminum  oxide.  The  0-rings  provide  a high-compliance 
suspension  to  the  inner  assembly  to  mechanically  isolate  and  seal  it 
inside  the  housing.  The  dielectric  housing  electrically  insulates  the 
pistons  which  act  as  electrodes  for  the  piezoceramic  tube.  The  housing 
also  acoustically  shields  the  inner  assembly  so  that  only  the  outside 
faces  of  the  pistons  are  exposed  to  sound  pressure.  The  sensor  element 
in  this  hydrophone  is  mounted  within  an  elastomer  boot  filled  with  cas- 
tor oil. 


HIGH  MASS  PISTON 
O'  RING 


CERAMIC  tubular 
HOUSING 


PIEZOCERAMIC  TUBE 

‘O'  RING 

LOW  MASS  PISTON 


Fig.  1.  Sensor  element  cross  section 

2 3 4 

Others  * * reported  on  hydrophone  sensors  of  this  configuration 

but  did  not  treat  the  directivity  aspect  of  this  design  in  the  vertical 

plane. 


Eleatromeahaniaal  Cirauit  Analysis 

The  physical  parameters  of  the  sensor  configuration  described  are 
shown  in  Figure  2.  When  the  piston  surfaces  are  exposed  to  a sound 
field,  the  inner  assembly  vibrates  about  a nodal  plane  perpendicular 
to  the  piezoceramic  tube  axis. 


2 

A.  A.  Aman'eva,  Ceramic  Acoustic  Detectors  (English  translation  pub- 
lished by  Consultants  Bureau  Enterprises,  Inc.,  New  York,  N.  Y.,  1956), 
Chapter  4,  pp.  80-83. 

^H.  A.  J.  Rijnja,  Acustica  30,  20-29  (1974). 

*H.  A.  J.  Rujnja,  Acustica  33,  1-9,  (1975). 
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Fig  or  element  physical  parameters:  nr]  is  the 

sm a'  mass,  m2  the  larger  piston  mass,  and  ap  the 

pi  us;  for  the  piezoceramic  tube,  mc  is  its  mass, 

lc  ts  length,  ai  its  inside  radius,  and  a2  its  outside 
radius;  l]  is  the  distance  from  the  lower  mass  piston  to 
the  nodal  plane,  I2  the  distance  from  the  higher  mass  piston 
to  the  nodal  plane,  and  ls  the  overall  sensor  length. 


The  electroacoustic  equivalent  circuit  is  shown  in  Figure  3,  where 
the  circuit  parameters  are  identified,  and  is  based  on  the  following 
assumptions: 


. That  it  operate  below  electroacoustic  resonance  where  the  sensor 
sensitivity  is  constant  with  frequency; 

. that  the  acoustic  radiation  and  mechanical  resistances  are  negli- 
gible; 

. that  it  operate  at  frequencies  below  the  mechanical  resonance 
of  the  pistons;  and 

. that  for  the  piezoceramic  tube 

$ = (s33ac)/(^33^c) » 

Cml  = sbD3*i/Sc  and  Cm2  = s3D3*2/Sc  (1) 


in  m/N,  and  Cg  = £(,5  /fc  in  ^ * where 

D 2 

s33  = open  circuit  elastic  compliance  in  m/N, 

d33  = piezoelectric  constant  in  m/V, 

Sc  = *(a22  - ax2)  in  m2,  and 

Ek  = blocked  permitivity  in  F/m. 

. Also,  if  we  assume  that  the  mechanical  resonance  frequency  of  the 
two  mechanical  branches  in  the  circuit  of  Figure  3 are  equal,  then 

"0  - ■ (‘/Wdz’1*-  '2> 

If 


mdl  = mi  + (mc/3)  (^Ag). 
md2  = m2  + (mc/3)  (i2/*c). 


and 


h + l2  = 


‘c* 


(3) 

(4) 


(5) 


substituting  Eqs.  (1),  (3),  (4),  and  (5)  into  Eq.  (2)  results  in 

mr  + mr/3 

m 1 « = m..  + ( — * — ) (— — — 7 — * — pr) 

dl  1 v 3 ' vm^  + m2  + 2m(./3/ 


(6) 


and 


md2  = + (-I-)  ( 


+ mc/3 


^ + m2  2mc/3 


-) 


(7) 
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k. 


* I 


Fig.  3.  Sensor  electroacoustic  equivalent  circuit:  Ce  and 
Re  are  the  blocked  electrical  capacitance  and  resistance,  $ 
is  the  voltage/ force  transduction  ratio,  Rmi  and  Rm2  and 
Cmi  and  Cm2  are  the  mechanical  resistances  and  compliances 
of  £1  and  respectively,  m^-j  and  mj2  are  the  dynamic 
masses  on  either  side  of  the  nodal.. plane,  F]  and  F2  are 
force  generators,  and  E0  is  the  open-circuit  sensor  voltage. 


Figure  4 shows  the  electroacoustic  equivalent  circuit  transformed  to 
the  equivalent  electrical  circuit  and  simplified  by  the  foregoing  assump- 
tions. A network  analysis  by  solving  two  simultaneous  loop  equations  for 
Eo  as  a function  of  radian  frequency  w gives 


Eo  = <j)F 


1 


Cml(1-U)2Cm2md2)  + Cm2  (F2/fV  ^1_aj2Cmlmdl) 


Cml  + Cm2  + 4>kCe“u,^('mlCm2^mdl  + md2)  -♦z“zcetcm2md2 
+ cmlmdl"“2CmlCm2mdlmd2) 


(8) 


where  u>  = 2-nf  in  rad/sec  and  f = frequency  in  Hz. 


— — 


J 


tk HtftrVi. 


Fig.  4.  Circuit  of  Fig.  3 transformed  and  simplified  to  give 
equivalent  electrical  circuit.  Parameters  same  as  Fig.  3 


Theoretical  Directivity 

This  hydrophone  configuration  is  omnidirectional  in  the  horizontal 
plane  due  to  symmetry.  and  F^  represent  the  amplitude  of  the  pressure 

forces  acting  over  the  areas  of  pistons  1 and  2 caused  by  a plane  acoust'c 
wave  traveling  at  angle  3 as  shown  in  Figure  2.  Appendix  A shows  the 
derivation  of  the  equations  for  these  forces.  Substitution  of  Eqs.  (All) 
and  (A17)  of  Appendix  A into  Eq.  (8)  gives  the  open  circuit  voltage  sen- 
sitivity of  the  sensor  element  for  small  ka  in  the  vertical  plane  (x-z): 
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Eq.  (9)  reduces,  for  small  values  of  ka„  and  la.  to  the  familiar 
form  p s 


Eo 

P 


<>C  S 
* m c 


The  theoretical  directivity  of  this  sensor  configuration  in  the 
vertical  plane  can  be  derived  by  taking  the  ratio  of  Eq.  (9)  evaluated 
at  angle  e to  that  at  e = 0°.  For  0 < e < 90°  one  obtains 


Dm  + cos  (ks.$  cos  e)  2J^  (kap  sin  e) 
D_  + cos  {kO  kap  sin  e 


where  M is  the  sensitivity  ratio  for  the  space  below  the-  x-y  plane 
relative  to  that  on  the  -Z  axis.  For  the  space  above  this  plane, 

90  < 9 < 180°  , one  obtains 


Dm  cos  [k«.$  cos  (180°  -e)]  + 1 2Jj  [ka  sin  (180°  -e)] 


Dm  + cos 
m 


kap  sin 


From  Eq.  (12)  it  can  be  seen  that  for  the  maximum  reduction  in  sensi- 
tivity above  the  x-y  plane  the  numerator  of  the  first  term  in  brackets 
of  Eq.  (12)  must  be  zero  at  6 * 180°.  Then  at  this  angle 

0m  cos  (k*s)  +1=0, 

and  solving  for  this  value  of  k£.$ 

(ka$)0  = Arcos  (-1/Dm).  (1 

Examination  of  Eqs.  (11)  and  (12)  shows  that  tr.e  left  term  in 
brackets  is  a function  only  of  Dm,  kz  , and  6.  The  right  term  in 

brackets  is  a function  only  of  ka  and  0.  Note  that  this  term  is 

identical  to  the  directivity  temrfor  a plane  piston.  Eq.  (11)  shows 

that  for  large  values  of  D and  small  values  of  ka„  both  terms  in 
3 m p 

brackets  approach  unity  and  the  sensor  element  will  approach  the  omni- 

directive  requirement  in  the  space  below  the  x-y  plane  (Figure  2). 

Figures  5 and  6 are  computer  generated  plots  using  Eqs.  (11)  and 
(12)  in  d8,  showing  the  effect  of  Dm  and  kap  on  the  directivity  from 

0°  < 9 < 180°  (x-z  plane  of  Figure  2).  Figure  5 shows  the  directivity 
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for  values  for  Dm  of  1 and  2.55  where  kap  = .787.  For  this  plot  (k«.s)g 
is  calculated  from  Eq.  (13)  for  each  Dm.  Figure  6 shows  the  directivity 
for  Dm  = 2.55  and  values  of  kap  of  0.1  and  1.0.  This  plot  shows  that 
there  is  very  little  effect  of  kap  for  values  up  1.0.  When  kap  becomes 
much  greater  than  1.0,  Eqs.  (11)  and  (12)  become  invalid  because  of 
assumption  that  A>>ap,  for  these  equations  (see  Appendix  A).  The  true 
cardioid  directivity  pattern  occurs  only  at  (k&s)0  as  can  be  seen  from 
these  plots. 


9 = 180* 


Fig.  5.  X-Z  directivity  patterns  of  sensor  element  where 
0°  < 9 < i80°  for  Dm  = 1 (dashed  line)  and  Dm  = 2.55 
(soTid  Tine)  calculated  from  Eqs.  (11)  and  (12); 
k£s  = (k£s)o  calculated  from  Eq.  (13)  and  kap  = 0.787. 
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Fig.  6.  X-Z  directivity  patterns  of  sensor  element  where 
0°  <_  4>  <_  180°  for  kap  = 0.1  (solid  line)  and  kap  = 1.0 
(dashed  line)  calculated  from  Eqs.  (11)  and  (127; 

Dm  = 2.55  and  k£s  = (k£s)o  calculated  from  Eq.  (13). 

These  directivity  patterns  are  shown  only  for  0°  < e < 180°  because 
for  180°  < 9 < 360°  the  pattern  is  the  mirror  image  of  that  shown.  This 
is  because  the  sensor  is  symmetrical  in  the  x-z  plane. 

Consideration  of  Eq.  (13)  shows  that  ( kx-s )q  can  have  an  infinite 
number  of  values  but  only  the  range  of  tt  < (k£s)g  < */2  corresponding 
to  1 < Dm  <»  are  of  practical  interest.  Eq.  (12)  shows  that  the  maximum 
reduction  in  sensitivity  occurs  at  e = 180°  for  any  value  of  Dm-  Also 
• the  sensitivity  at  9 = 180°  increases  as  k^  changes  in  either  direction 

from  (k*s)Q  within  the  above  range.  If  e = 180°  and  for  small  values 
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2J,  (ka  ) 

of  kap  where  y — = 1,  the  values  of  kJ.$  on  either  side  of 

(kzs)g  can  be  foundp  where  < 0.5  (-6  dB).  For  these  conditions  Eq. 
(12)  becomes 

D cos  kx.  + 1 

Me  = D + cos  kx  3 ±0.5  . 
m s 

Solving  for  kx„  for  the  two  values  of  Ma 

(kz$)1  = Arcos  [ (Dm-2)/(2  Dm  +1)]  and  (14) 

(kz$)2  = Arcos  [-(Dm  +2)/ (2  Dm  +1)].  (15) 

Figure  7 shows  Eqs.  (13),  (14),  and  (15)  plotted  for  comparison.  The 
plot  of  (kx.$)g  shows  that  it  rapidly  decreases  from  w radians  (x.$  = 

x/2)  as  Dm  increases  from  1.0  and  approaches  ir/2  radians  (x_s  * x/4) 

asymptotically.  The  plots  of  (kx^j  and  ( k^s )£  show  the  same  trend 

except  they  approach  tt/3  radians  (t  = X/6)  and  2it/3  radians  (x$  = 

x/3)  respectively  as  Dm  increases. 


Fig.  7.  Plots  of  Eqs.  (13),  (14),  and  (15).  The  area 
bounded  by  the  curves  for  (k£$)i  and  (k£s)2  contain  all 
values  for  k£s  and  Dm  where  M0  <_  0.5  (-6  dB)  at  e = 180°. 


Experimental  Direativity  Comparison 

The  validity  of  the  theoretical  directivity  developed  for  this 
sensor  configuration  was  verified  by  computer  generated  directivity 
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patterns  of  the  USRD  type  H56  hydrophone  using  Eqs.  (11)  and  (12).  Com- 
parison of  these  theoretical  plots  with  the  actual  directivity  measured 
at  the  USRD  Lake  Facility  is  shown  in  Figures  8a  thru  c for  three  fre- 
quencies. The  center  frequency  of  27  kHz  was  calculated  from  Eq.  (13) 
where  for  this  hydrophone  D = 2.46,  the  sensor  length  i = 17.45  run, 
and  ap  = 7.87  mm.  The  plots  show  good  agreement  be tween smeasu red  and 

calculated  plots.  The  worst  agreement  shown  was  for  a frequency  of  30 
kHz.  This  can  be  attributed  to  diffraction  due  to  proximity  of  the  hy- 
drophone preamplifier  housing  at  the  top  and  the  sensor  frame  ring  at  the 
bottom  of  the  pattern  and  the  fact  that  at  this  frequency  ka_  is  almost 
1.0.  p 


9 = ISO' 


Figs.  8a-c.  X-Z  directivity 
patterns  for  USRD  type  H56 
hydrophone  where  0°  < 0 < 180° 
for  3 frequencies  shown; “sol id 
lines  are  measured  patterns 
and  dashed  lines  are  patterns 
calculated  from  Eqs.(ll)  and 
(12). 


Fig.  8c 


After  verification  of  this  design  analysis  with  the  actual  direc- 
tivity of  the  H56  hydrophone,  the  USRD  type  H73  hydrophone  sensor  element 
was  developed  to  fulfill  the  requested  design  requirements.  Eleven  sensor 
elements  with  values  of  D from  1.8  to  10.4  were  evaluated.  After  com- 
parison of  the  calculatedmdirectivity  patterns  with  the  measured  patterns 
it  was  found  that  Eqs.  (11)  and  (12)  were  inaccurate  for  values  of  D 
above  2.55.  Figures  9a  thru  c show  the  type  H73  sensor  element  that 
gave  the  best  correlation  between  the  calculated  and  measured  direc- 
tivity patterns.  For  this  sensor  Dm  * 2.55,  ts  = 27.1  nm,  and  ap  * 11.1 

mm.  It  consists  of  an  8.51-mm  00  x 5.02-mm  ID  x 12.70-mm  long  PZT-5 
lead  zirconate-titanate  piezoelectric  tube.  The  bottom  piston  made  of 
type  6061-T6  aluminum  alloy  is  6.35-mm  thick  and  has  a mass  of  6.11  g. 

The  top  piston,  made  of  type  304  stainless  steel,  is  8.05-mm  thick  and 
has  a mass  of  23.05  g. 
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The  maximum  reduction  in  sensitivity  above  the  x-y  plane  occurred 
at  18  kHz.  The  calculated  frequency  for  this  to  occur,  from  Eq.  (13), 
was  17  kHz  where  ( ki-s )q  3 1.97.  The  agreement  between  calculated  and 

measured  directivity  patterns  is  not  as  good  as  that  for  the  type  H56 
hydrophone.  The  open-circuit  free-field  voltage  sensitivity  of  the 
sensor  in  this  frequency  range  was  -182  dB  re  1 V/u Pa. 

Conclusions  And  Recommendations  . 

Based  on  the  analysis  verified  in  this  report  a hydrophone  sensor, 
such  as  the  USRD  type  H73  hydrophone,  can  be  configured  that  is  omni- 
directional in  the  space  below  the  sensor's  horizontal  plane  and  has 
reduced  sensitivity  in  the  space  above  this  plane.  Once  the  parameters 
are  established  for  this  configuration,  its  characteristic  directivity 
pattern  will  approach  a true  cardioid  shape  at  only  one  frequency.  The 
desired  reduction  in  sensitivity  above  the  sensor's  horizontal  plane 
occurs  in  a limited  frequency  range. 

The  procedure  in  designing  a sensor  element  of  this  configuration 

knowing  its  desired  directivity,  its  frequency  range,  and  its  dynamic 

mass  ratio  D„  is  as  follows: 
m 

1)  Using  Eq.  (13)  or  Figure  7 calculate  (k£  )Q  where  the  theoretical 
cardioid  directivity  pattern  occurs.  s 

2)  Calculate  the  sensor  length  i from  (kt  )Q  at  the  center  of  the 
desired  frequency  range. 

3)  Using  Eqs.  (14)  and  (15)  or  Figure  7 calculate  (k£$)j  and  (kJi^. 

4)  Using  Eqs.  (11)  and  (12)  calculate  the  theoretical  directivity 
patterns  at  other  frequencies. 


This  procedure  can  easily  be  adapted  to  be  carried  out  on  a digital  com- 
puter with  the  directivity  patterns  drawn  by  an  x-y  plotter.  This  anal- 
ysis is  valid  for  values  of  ka  <1.0  (a  < .16a)  and  values  of  dynamic 
mass  ratio  Dm  < 2.55.  p p 

The  limitation  of  this  analysis  to  the  stated  values  of  D should 
be  investigated  to  determine  the  cause.  This  limitation  is  prBbably 
caused  by  the  assumption  to  neglect  the  acoustic  impedance. 
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APPENDIX:  WAVE  ANALYSIS 


The  forces  and  in  Eq.  (8)  represent  the  force  amplitudes 

caused  by  the  sound  pressure  acting  on  pistons  1 and  2 (represented  by 
masses  m^  and  in  Figure  2).  The  force  on  piston  I for  x»ap  is 

given  by  the  equation 


Fj  (e,t)  ■ 4 pj  (x,9,t)  ydx  (Al] 

where  p^  (x.e.t)  is  the  sound  pressure  distribution  on  the  piston  and 

ydx  is  the  differential  area  over  which  p,  (x,e,t)  acts  at  any  coor- 
dinate x.  1 

The  plane  wave  equation  on  the  Z*  axis  is 


p (z',t)  = pp  (cos  ait) (cos  kz')  (A2) 

where  pp  is  the  sound  pressure  amplitude,  u and  k are  as  defined  before, 

and  z'  is  the  distance  from  the  origin.  For  any  coordinate  x on  piston 
1 the  corresponding  coordinate  along  the  Z'  axis  (Figure  1A)  is 

z'  = x sin  e.  (A3) 


Fig.  1A.  Coordinate  system  for  sensor  looking  at  piston  1. 
Substituting  Eq.  (A3)  in  Eq.  (A2)  gives 

Pj  (x.e.t)  s p (cos  wt)  [cos  (k  x sin  e)]. 


i II  Hi  i ' 


A-l 


Looking  at  Figure  1A,  it  can  be  seen  that 
y - (ap2  - X2)15 


and  therefore 


ydx  = (ap2  - x2)15  dx. 


Also  note  from  Figure  1A 

x = a„  sin  6 and  dx  * a„  cos  8d  6. 
P P 

Substituting  Eq.  (A6)  into  Eqs.  (A4)  and  (A5)  gives 


p.|(x,e,t)  = p^B.e.t)  * pp(cos  wt)cos[(kapsin  e)(sin  b)],  and 

2 2 

ydx  = ap  (cos  b)  dB. 


Substituting  Eq.  (A7)  into  Eq.  (Al)  gives  the  final  transformation  to 
polar  form  ^2 

2 2 

F,(0,t)  = 4p  a cos  ait  (cos  b)  cos[(ka  sin  e)(sin  s)]  dB.  (A8) 
1 P P , ^ P 


From  Reference  5 the  following  definite  integral  expression  for  a 
Bessel  Function  of  order  n where  n>-*s  is  given  as 

J (u)  = (7)°  (cos  a)2n  cos(u  sin  a) da.  (A9) 

n /T-r(n+5s)  * Jo 

If  n = 1,  a = 6 and  u = ka  sin  0 and  knowing  that  r(n  + *5)  = */2  Eq. 
(A9)  becomes  p 


J^kapSine) 


2kapsine  .n/2 


(cos  B)2  cos[(ka_sin  0)  (sin  8)]dB.  (A10) 
Jo  p 


Substituting  the  term  on  the  left  side  for  the  terms  in  the  right  side 
of  Eq.  (A10)  contained  in  Eq.  (A8)  gives 


2J.  (ka  sin  0) 
Fl(e,t)  * pp  sp  cos  „t  - kapi?n- 


T.  M.  McRoberts,  Spherical  Harmonics  (Dover  Publications,  Inc.,  New  York, 


N.Y.,  1947)  Chap.  14,  p'p  255^7. 


where  Sp  = iap  , the  piston  area. 


Finally,  the  amplitude  of  F^e.t)  is 


Fi  • pp  sp 


2J1  (kap  sin  8) 
ka_  sin  e 


(All) 


The  force  on  piston  2,  for  x»a  is  given  by 

P 


F2(6,t)  = 4 qP  p2(x,9,t)ydx 


(A12) 


where  p2(x,e,t)  is  the  sound  pressure  distribution  on  this  piston 
and  ydx^is  differential  area  over  which  p2(x,e,t)  acts  at  any  coor- 
dinate x.  For  any  coordinate  x in  piston  2 the  corresponding  coordi- 
nate along  the  Z'  axis  (Figure  1A)  is 


z1  = i cos  9 + x sin  0. 


(A13) 


Substituting  Eq.  (A13)  in  Eq.  (A2)  gives 

p2(x,0,t)  = Pp(cos  ut)  cos[k(ts  cos  0 + x sin  0)].  . (A14) 

Transforming  Eq.  (A14)  and  by  substituting  for  x from  Eq.  (A6)  gives 

P2(x,9,t)  * p2(s,0,t)  = pp(cos  a)t)cos{k[£scos  0-ap(sin  0 ) ( si n 6)]}.(A15) 

Substituting  Eq.  (A15)  and  with  the  area  ydx  from  Eq.  (A7),  Eq.  (A12) 
becomes 


F2(0,t)  = 4ppap  cosuit 


r n /2  y 

(cos  e)  cos{kLficcos  0 + a (sin  0)(sin  s) ] }dB . 


This  equation,  using  a trigonometric  identity,  becomes 

i-tt/2 


F2(e,t)  * 4ppap  (coswt) [cos (k)is cose)] 


(cose)  cos[ka  ( sine) (sine) ]d6 
0 p 


2 t^/2  n 

-4ppap  (cosu)t)[sin(kK.scose)]  (cose)  si n[kap(sine)( sine) ]de. (A16) 


The  second  term  of  Eq.  (A16)  when  integrated  is  zero.  Finally,  the  first 


term  in  Eq.  (A16),  after  rearranging  to  the  form  of  Eq.  (A10) , gives 


